Anisotropy of the Mobility of Pentacene from Frustration 
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The bandstructure of pentacene is calculated using first-principles density functional theory. A 
large anisotropy of the hole and electron effective masses within the molecular planes is found. The 
band dispersion of the HOMO and the LUMO is analyzed with the help of a tight-binding fit. The 
anisotropy is shown to be intimately related to the herringbone structure. 
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I. INTRODUCTION 

Charge transport in organic molecular conductors has 
gained prominent interest as it is expected that in ultra- 
clean materials band transport with high electronic mo- 
bilities can be obtained. Recent progress in synthetic 
procedures have minimized mobility reducing factors, 
such as disorder, Frenkel defects, and Schottky defects. 
Thus, mobilities up to ~ 1 cm 2 /Vs have been ob- 
tained in thin film devices and in single crystals at room 
temperature. 1-4 

Pentacene has gained prominent interest because of a 
number of favorable materials properties. Large single 
crystals can be easily grown by vapor transport, 5 and 
the material can be handled in air. The reactivity of pen- 
tacene to impurity levels of oxygen and hydrogen is often 
underestimated, but can be minimized using appropriate 
synthetic conditions. 6 Despite the presence of at least 
four different polymorphs, 7 there are no structural phase 
transitions observed for the most stable polymorph, the 
single crystal polymorph. 8-10 The layered crystal struc- 
ture indicates a highly two-dimensional charge transport. 

In this paper we use first-principles density functional 
theory (DFT) to calculate the bandstructure of pen- 
tacene. We find a very large anisotropy of the effective 
mass within the layers, resulting from a frustration of the 
electron transfer along the crystallographic b-axis. This 
frustration is rationalized using a tight-binding fit of the 
band structure, and is intimately related to the herring- 
bone geometry of the crystal structure of pentacene. Our 
results imply that the mobility is strongly directionally 
dependent within the layers. 



II. STRUCTURE 

The crystal structure of bulk pentacene, depicted in 
Fig. 1, has been reported by several groups. 8-10 It is 
a layered structure, triclinic, with spacegroup PI. The 
molecular layers consist of pentacene molecules arranged 



in a herringbone pattern. We chose the unit vectors fol- 
lowing Ref. 9, with a and b in the molecular plane and 
c inter- relating the molecular planes. The unit cell con- 
tains two inequivalent molecules. The molecules are not 
perpendicular to the molecular planes, but are at finite 
angles. This slant can be understood as a relative shift 
of the molecules. Starting from one molecule, moving 
along (a + b)/2 one arrives at an another (inequivalent) 
molecule, that is shifted along its longest molecular axis 
(LMA) by 2.13 A, i.e. approximately the length of one 
ring in the pentacene molecule. Going along (a— b)/2 the 
shift is a mere 0.12 A, i.e. it is negligible. This "staircase" 
pattern is indicated in Fig 1(b), which shows a projection 
of a molecular layer along the LMA. 



III. TECHNICAL DETAILS 

The bandstructure has been calculated using den- 
sity functional theory in the generalized gradient ap- 
proximation (GGA). 11 The ab initio total-energy and 
molecular-dynamics program VASP (Vienna Ab initio 
Simulation Program) was used. 12-15 Electron- ion inter- 
actions were described using the projector augmented 
wave method. 16 ' 17 The kinetic energy cutoff on the wave 
function expansion was 500 eV. 

The self-consistent calculations were carried out us- 
ing the experimental positions and cell (determination 
at 90 K from Ref. 9) with a 4x4x2 k-point mesh. 18 Con- 
vergence was tested with a 6x6x4 mesh and found to be 
sufficient. 

The bandstructure of the single molecular layer was 
calculated also from a 3D crystal, but with an extra vac- 
uum region of 6 A thickness inserted between the lay- 
ers. Thus interactions between the layers were avoided. 
Moreover, the c-axis was chosen perpendicular to the ab 
plane, i.e. the lattice was transformed to that of a mon- 
oclinic crystal. 
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IV. RESULTS 

Fig. 2 displays the calculated density of states (DOS) 
of pentacenc. All bands below the Fermi level are com- 
pletely filled. The unoccupied bands are separated from 
the valence bands by a band gap of 0.7 eV. This is consis- 
tent with experiment, as undoped pentacene is an insu- 
lator. The experimental band gap is 2.2 eV. 19 However, 
an underestimation of the band gap is not unusual for 
GGA/DFT and, in general, doesn't much affect the other 
features of the electronic structure. The DOS peak of the 
HOMO (just below eV) and of the LUMO (at 1 eV) 
are both well separated from the other peaks. Each of 
them doesn't mix with bands derived from other molec- 
ular states. The calculated HOMO and LUMO band 
widths are ~ 600 and ~ 700 meV respectively. 

Another, recent theoretical calculation gives values of 
608 and 588 meV respectively. 20 This compares well with 
our results, although we find a bit wider LUMO band. 

Both the DOS of the single layer and of the complete 
3D crystal are plotted. On the scale of the plot, differ- 
ences are only apparent above 3 eV. The electronic struc- 
ture of the single layer should be a good first approxima- 
tion to that of the complete crystal. Its bandstructure is 
only 2D, which greatly facilitates the interpretation. 

A. Bandstructure of a Single Layer 

The band structure near the Fermi level of a molecu- 
lar layer is depicted in Fig. 3. Both the HOMO and the 
LUMO complex consist of two bands, as there are two 
molecules in the unit cell. For both, the maximum and 
minimum are at M [(5,5)]. Therefore, upon doping (both 
with electrons and with holes) the charge carriers will be 
created in the vicinity of this point. Their effective mass 
is strongly directionally dependent. In particular, note 
the flatness of the bands near the extrema for the direc- 
tion M-Y. These observations will be further elaborated 
upon below. 

Since the band complexes are well separated from the 
other bands, we can attempt a tight-binding (TB) fit, 
for both the HOMO and the LUMO separately. The 
fit parameters are the transfer integrals t between the 
molecular orbitals. Some of these are shown in Fig. 4. 
As we have two molecules, a 2 x 2 matrix will have to be 
diagonalized: 

"e(k)-e-Ti T 2 ] 

T 2 £(k)-7ij [i) 

Here T\ and T 2 denote the sums of the diagonal and off- 
diagonal transfer integrals, multiplied with the appropri- 
ate phase factors, respectively. e(k) are the eigenvalues 
and e accounts for the fact the the molecules are inequiv- 
alent. e(k) is required to follow the DFT result of Fig. 3 
as close as possible. 



As is evident from the figure, the fit follows the DFT 
result very well. All transfer integrals are listed in Ta- 
ble I. By far the most important ones are t ai t( a +b)/2 
and i( a -b)/2- For the HOMO (LUMO) the values are: 
t a = 31 (-41) meV, t( a +b)/2 = -56 (-90) meV and 
£(a-b)/2 = 91 (90) meV. All the others are very small. 

The effective mass is easily calculated with the help 
of the analytic expression for the fit. In Fig. 5 the hole 
effective mass m*, calculated from: 

d 2 e(k) = m 
Ok 2 to* [ ' 

is plotted for all directions in the plane. Note the large 
anisotropy. The largest value, 58 m c , occurs for the a 
direction. 

Below we'll further discuss the band structure near the 
M point (i.e. the point where the valence band maxi- 
mum and conduction band minimum occur) and explain 
the mechanism behind the large anisotropy. The molec- 
ular orbitals, HOMO and LUMO, are plotted in Fig. 6. 
They are both 7r-systems, having a nodal plane coincid- 
ing with the plane of the molecule. Along the LMA the 
HOMO alternatingly has maxima and minima that are 
slightly more apart than the length of one ring-unit of the 
molecule. This pattern leads to 4 nodal surfaces "per- 
pendicular" to the LMA. Additionally there's one more 
nodal plane containing the LMA. 

The LUMO doesn't have such a nodal plane. It has two 
more nodal surfaces "perpendicular" to its LMA. Again 
the wave function extrema are approximately one ring- 
unit apart. 

Fig. 7 shows a schematic representation, in a plane 
perpendicular to the LMA, of the HOMO-derived crystal 
orbitals at the M point. At M the phase of the wave 
function changes by 180°, both for a translation over a 
and b. However, in Fig. 7 all 4 molecules at the corners 
have the same phase. This is because the staircase (slant) 
in the layer gives rise to a shift of the molecule at each 
step. This results (see the preceding paragraph) in an 
additional phase change of approximately 180°. 

Depending on the phase of the central molecule, either 
the low or high band at M is selected. In Fig. 7(b) all 
three important transfer integrals give a bonding con- 
tribution. Indeed, this is the most bonding eigenstate 
within the band complex and thus corresponds to the 
minimum at M. Fig. 7(a) depicts the most antibonding 
situation (the top of the valence band): The transfers 
(a + b)/2 and (a — b)/2 arc maximally anti-bonding. 
Transfer a is still bonding. Due to frustration no solu- 
tion with anti-bonding overlap for all three transfers is 
possible. 

The frustration can also explain the large anisotropy in 
the effective mass. Moving along a (« a*) two conduction 
channels are possible: One directly via t a , the other via a 
"double hop" [f( a +b)/2 and i( a -b)/2]- Moving away from 
M, the double hop will force the bands in a downward 
curvature, which is counteracted by the direct channel. 
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This explains why the band becomes flat. By contrast, 
moving along b* in reciprocal space, the channel along a 
drops out and a cosine-shaped band results, indicative of 
alternate (a + b)/2 and (a — b)/2 hops. 

The LUMO bands mimic the HOMO bands in a very 
special way: They appear as the "mirror" image of the 
HOMO bands on the energy axis (Fig. 3). This pecu- 
liar feature is due to the LUMO's slightly different nodal 
structure (it lacks the second mirror plane containing the 
LMA). Therefore the top point at M has all three trans- 
fers antibonding, whereas frustration now occurs for the 
most bonding point at the bottom of the conduction band 
(where the a "hop" is still antibonding). Thus the large 
anisotropy of the hole effective mass is "mirrored" in the 
electron effective mass. The analogy between the HOMO 
and LUMO band complex is not complete: E.g., the 
avoided crossing along M-r only occurs for the HOMO 
and the total width of the LUMO complex is slightly 
larger than that of the HOMO complex. 

B. Bandstructure of the 3D Crystal 

The band structure of the full 3D crystal (Fig. 8) 
strongly resembles that of the 2D layer. The shape of 
the bands for the most dispersive, a* and b* directions 
is essentially the same. Some differences are apparent 
though: (a) the top of the valence band is not at M 
[(i,i,0)], but at (0.375, 0.5, 0.075). Moreover, as is evi- 
dent from the top panel of Fig. 8, the band doesn't have 
a parabolic shape. The deviations are small however. 

The conduction band minimum now occurs near H. (b) 
A small band dispersion occurs along c* . 

Fig. 9 depicts the hole effective mass, directly calcu- 
lated from the band structure using finite differences. We 
cannot rely on a TB fit, as it cannot capture the bump at 
(0.375, 0.5, 0.075). For directions close to a* the effective 
mass was not calculated, due to the locally non-parabolic 
nature of the band. For the other directions the picture 
is very similar to the 2D analogon. For c* we obtain: 
to*. = 5.2 m e . 

Our results may be compared with a recent extended 
Hiickcl (EH) band structure calculation on the pentacene 
single crystal. 21 The agreement between the EH and our 
first-principles results is quite good. The EH study re- 
ports the band structure for a limited number of lines in 
the BZ (r-X, T-Y, T-Z) and therefore doesn't consider 
the special situation near M. 

V. DISCUSSION & CONCLUSIONS 

Band structure calculations for the pentacene single 
crystal reveal HOMO and LUMO bandwidths of ~ 600 
and 700 meV respectively. Consistently, (hole) effective 
masses as low as 1.7 m e (for m{,«) are found. Even to*. 
(5.2 m e ) is rather low, in spite of the apparent flatness 



of the bands for the c* direction. Within the molecular 
layers a competition between different conduction paths 
results in a large anisotropy of both the hole and the 
electron effective mass, the charge carriers being orders of 
magnitude heavier for the a* (than for the b* direction). 

Our calculations do not take into account the cou- 
pling with lattice degrees of freedom. At lower tem- 
peratures the electron-phonon coupling may renormal- 
ize the charge-carrier's bandwidth. At higher tempera- 
tures the bandwidth renormalization and mobility reduc- 
tion may differ for different directions, thus affecting the 
anisotropy. For anthracene, recent first-principles calcu- 
lations indeed show a directionally dependent bandwidth 
renormalization. 22 However, as it is so pronounced, we 
expect much of the effective mass anisotropy to survive, 
also at higher temperatures, as long as coherent band- like 
transport is dominant. 
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TABLE I. The following table lists the parameters for the 
single layer TB fits. 





HOMO 


LUMO 


e (eV) 


0.042 


0.047 


t a 


0.061 


-0.041 


h 


-0.004 


-0.008 




-0.001 


0.000 


t a + b 


0.003 


-0.002 


t a -b 


0.000 


-0.002 


t(a + b)/2 


-0.056 


-0.090 


t(a-b)/2 


0.091 


0.090 


*3(a+6)/2 


-0.001 


0.000 


*3(a-6)/2 


0.001 


0.000 


*(3a+b)/2 


0.000 


-0.001 


t{3a-b)/2 


0.001 


0.001 


tae 


0.000 


0.003 




FIG. 1. The crystal structure of pentacane. (a) Stacked 
layers of pentacene molecules, viewed along the [110]-axis. A 
unit cell is also indicated, (b) Projection of the pentacene 
crystal structure and its unit cell vectors a and b on a plane 
perpendicular to the long molecular axes (LMA). The angle 
between the molecules is indicated. The herringbone arrange- 
ment is evident. 
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FIG. 2. Density of states in states per eV per cell as func- 
tion of the energy. The energy is chosen zero at the top of the 
valence band. Solid line: single crystal, dotted line: cell with 
the molecules 6 A apart. Both curves are calculated using a 
4x4x2 mesh. IDOS is the integrated DOS in states(electrons) 
per cell. 




FIG. 4. The pentacene afo-plane, projected along the long 
molecular axis. The directions of the different overlap inte- 
grals are indicated. 




r X M Y r M r FIG. 5. The reciprocal effective mass as function of the 

crystallographic direction. Calculations are based on the 
tight-binding fit of the band structure, which is calculated for 
FIG. 3. The band structure of one layer of pentacene one layer of pentacene molecules. The pentacene molecules 

molecules, with a layer spacing of 6 A. The energy is cho- are shown as reference to the crystallographic directions. 

sen to be zero at the top of the valence band. The energies 

are plotted along lines in the first Brillouin zone, connect- 
ing the points F = (0,0), X = (a*/2,0), M = (a*/2,6*/2), 

and Y = (0, 6*/2). From T to M is in the direction of 

l/2(a* + b*), from M to Y is from l/2(a* - b*) to Y. The 

solid line is the tight binding fit. 
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HOMO 




LUMO 




FIG. 6. Calculated electron densities of the HOMO and the 
LUMO in a plane above the molecular plane. The large and 
small black circles indicate the position of the C-atoms and 
H-atoms, respectively. The lines are on a logarithmic scale. 
The lines lie at 0.0001 • 10 n/3 elec/A 3 , n > 1. The outermost 
lines have the same electron density. 
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FIG. 8. The band structure of single crystalline pen- 
tacene The energy is chosen to zero at the top of the va- 
lence band. The energies are plotted along lines in the 
first Brillouin zone, connecting the points T = (0,0,0), 
X = (a*/2,0,0), M = (072,672,0), Y = (0,672,0), 
Z = (0,0, c72), TV = (0,672, c72), H = (a* /2,b* /2,c* /2) 
and O — (a*/2, 0, c*/2). From r to M is in the direction of 
l/2(a* + b*), from M to Y is from l/2(a* - b*) to T. The 
solid line is the tight binding fit. The top graph shows an 
enlargement of the top of the valence band. 
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FIG. 7. Crystal orbitals of the HOMO complex at M 
[(5,5)], viewed along the LMA, in a plane perpendicular to 
the LMA. The HOMO wave functions are schematically de- 
picted with the + and - signs to indicate the phase, (a) The 
most antibonding situation, as is the case at the top of the 
valence band, (b) All orbital overlaps are bonding, forming 
the bottom of the valence band. The dashed lines indicate 
the shifts between the molecules. 



fa 




270 

FIG. 9. The reciprocal effective mass as function of the 
crystallographic direction. Calculations are based on the 
band structure of single crystalline pentacene. The pentacene 
molecules are shown as reference to the crystallographic di- 
rections. 
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